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Abstract 



A Condorcet voting scheme chooses a winning candidate as one 
who defeats all others in pairwise majority rule. We provide a review 
which includes the rigorous mathematical treatment for calculating 
the limiting probability of a Condorcet winner for any number of can- 
didates and value of n odd or even and with arbitrary rank order 
probabilities, when the voters are independent. We provide a com- 
pact and complete Table for the limiting probability of a Condorcet 
winner with three candidates and arbitrary rank order probabilities. 
We present a simple proof of a result of May to show the limiting 
probability of a Condorcet winner tends to zero as the number of can- 
didates tends to infinity. We show for the first time that the limiting 
probability of a Condorcet winner for any given number of candidates 
m is monotone decreasing in m for the equally likely case. This, in 
turn, settles the conjectures of Kelly and Buckley and Westen for the 
case n — > oo. We prove the validity of Gillett's conjecture on the 
minimum value of the probability of a Condorcet winner for m = 3 
and any n. We generalize this result for any m and n and obtain 
the minimum solution and the minimum probability of a Condorcet 
winner. 
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1 Introduction 



Different voting procedures (see for example [22] Page 147) for electing a 
leader or deciding a candidate or a judgment have been proposed in the 
literature. In all voting procedures, we assume that the voter rank orders 
the candidate or decision from highest to lowest. One of the most used 
schemes is a majority rule, in which a candidate or decision is ranked first by 
more than half the number of voters. An example of the majority rule is the 
U.S. Supreme Court decisions. Another selection strategy is the plurality 
rule in which a candidate is ranked first by the largest number of voters. 
An example of the plurality rule is the U.S. presidential election. In the 
Condorcet scheme, a candidate is chosen who defeats all others in pairwise 
majority rule. An example of the Condorcet scheme is the selection of job 
candidates or proposals to be funded in academic departments and by funding 
agencies or for determining winners in tournaments. 

While the Condorcet scheme is a transitive in individual voter's choices, 
there are many instances in which a majority winner may not exist [23] |19j . 
For example, let there be 3 voters (Vi, V2, V3) and 3 candidates (C\, C%, C3). 
V\ rank orders C1C2C3, V2 rank orders C2C3C1 and finally V3 rank orders 
C3C1C2. In this voting example, C\ is preferred over C2 two out of three 
times. Similarly C2 is preferred over C3 and C3 is preferred over C\ two out 
of three times. Hence there is no one with majority preferences over all pairs. 
So according to the Condorcet scheme, there is no winner in this case. 

In past research work, it has been assumed that a probability distribution 
is given to all the rank orders. It is also assumed in most studies that all vot- 
ers vote independently. Previous researchers have calculated the probability 
of a Condorcet winner using analytical and simulation models [T7| by mak- 
ing simplifying assumptions on the probability distributions. Asymptotic 
expressions for the limiting probability of a Condorcet winner (i.e., when the 
number of voters — > 00) have been obtained J2] [H] [22] P^j - A survey 
of previous research in this area may be found in |lUj . 

The focus of the present paper is on the limiting probability of a Con- 
dorcet winner when the number of voters, n, tend to infinity for a given 
number m of candidates and arbitrary probability distribution on the ml 
rankings. We present a review and some new results in this part of the study. 
While researchers in this area, notably for example, Niemi and Weisberg [2H] , 
Garman and Kamien [U] have noted the relationship of the limit with certain 
multivariate normal probabilities, there doesn't exist a statistically rigorous 
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proof of this result in order to justify the evaluation of the limit for various 
cases. Further it is often assumed that the number of candidates is odd, see 
Niemi and Weisberg We close this gap with this paper by providing a 
rigorous statistical proof for any number of candidates, the number of voters 
both odd and even, and arbitrary probability distributions, usually referred 
to as culture probabilities. In section 2, we show further that for any number 
of candidates and with arbitrary probability distributions on the preference 
ranking, the limiting probability depends on the calculation of positive or- 
thant probabilities of appropriate multinormal distributions. We also provide 
a compact and complete Table covering all possible cases for computing the 
limiting probability of a Condorcet winner when the number of candidates 
is 3 and with arbitrary probability distribution among the rank orders. We 
show that the limiting probability as n — > oo for any number of candidates m 
is monotone decreasing in m which partially validates Kelly ^H] and Buckley 
and Westen jl] conjectures. We also present a simple proof of May's result 
[T§] that the limiting probability of a Condorcet winner for equally likely 
case tends to zero as the number of candidates gets larger. We also prove 
Gillett's conjecture on the minimum probability of a Condorcet number 
for 3 candidates and any number of voters. We extend the result to any m 
and n and obtain the minimum value and solution. In Section 5, we treat 
the case of m = 4 candidates with arbitrary probabilities Pi and show that 
it is possible to have exact expressions for the limiting Condorcet winner 
probability for all possible scenarios. 

2 Condorcet Winner: Statistical Formulation 
and Results for General Case 

In this section, we present a statistical derivation of the limiting probability 
of a Condorcet winner when the voters are independent and the number of 
voters is large. Let 

m = Number of Candidates 
n = Number of Voters 

We have K = m\ preference rankings of m candidates. A voter will choose 
one of these rankings. We assume that the voters act independently. Let Pi 
be the probability that a voter prefers the rank order i, for i = 1, 2, ■ • • , K. 
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K 

Further, we know that Pi > and J^Pi = 1. Let iVj be the number of voters 
voting for the zth preference ranking, for % = 1, 2, • • • , K. We therefore have 

K 

/,Nj — n. Now, it is well known that (N\, N 2 , • ■ ■ , Nk) has a multinomial 

i=l 

distribution. 

nl * 

P(JVi = i = 1, • • • , JO = -p-p r II VT (1) 

K 

where the rij's are nonnegative integers with rij = n. 

i=l 

It is known that (see, for example, Page 318), for % = 1, • • • , K, j = 
1, • • • , K and i 7^ j. 



E[Ni] = n Pl , 
Var[Ni] = npi(l-pi 
Covariance(Ni, Nj) = —npiPj, 



Correlation(Ni, Nj) 



PiPj 



(l- ft )(l- Pi ) 



Denote the m candidates by Ci, C2, ■ ■ ■ , C m . Write CiPCj, i 7^ j if Cj is 
preferred to Cj by a voter. Define CiMCj, i 7^ j if Cj beats Cj by majority 
rule, i.e., the total number of voters preferring Cj to Cj is more than the 
total number of voters preferring Cj to Cj. This fact can be mathematically 
expressed as 

ai,(i,j) N i + a 2,{i,j)N 2 H h aK,(i,j) N K > 1 (2) 

where the a^j^'s are ± 1. a^jj) = 1, if in the /th preference ranking voted 
by the N voters, CiPCj. Similarly a^j) = — 1 if CjPCi is the preference 
ranking voted by iVj voters. Clearly ~ of the a^j^'s are 1 and y of the 
a «,(j,i)'s are equal to -1. 

Equation (J2J) implies that the difference between the number of voters 
preferring Cj to Cj and the number of voters preferring Cj to Cj is at least 1. 
We define, for example, Cj to be the Condorcet winner if and only if CiMCj 
for each j = 1, 2, • ■ ■ , m and j 7^ i. For example, the probability that C\ is 
the Condorcet winner is the probability of the joint event 

01,(1,2)^1 + 02,(1,2)^2 H V aK,(i,2)N K > 1 
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01,(1,3)^1 + 02,(1,3)^2 H h a K ,(i,3)N K > 1 

■ • ■ > 1 
... > 1 

Ol,(l,m) 02,(l,m) 

iV* > 1 (3) 

Condition (jSJ) lead to what is denned as 'strong' winner. Replacing all the 
l's by O's lead to the definition of a 'weak' winner. For small values n and 
m , this probability can be computed using the multinomial distribution (Q). 
We are interested in finding 

lim^oo P(One of the Cj is a Condorcet winner) for a given value of m. 

Several researchers have attempted to find this limit or an approximation 
to it. Gilbaud [H] gave an expression of this limit for m — 3, when all the 
Pi's are equal (Impartial Culture or IC model), though he did not explain 
the method. Garmen and Kamien j2] derived the expression for the limiting 
Condorcet winner probability for m = 3 and m = 4 for the Impartial Culture 
model without giving details. Jones et al [T7|, Bell j2] and Niemi et al 
[2Sj obtained approximations to it by extensive simulations or quadrature 
methods attributed to Ruben (25] for the Impartial Culture (IC) model i.e., 
when all the preference rankings are equal. 

Niemi and Weisberg j2H| and Garman and Kamien |JJ were the first to 
relate the limit or an approximation to it with the calculation of certain 
multivariate normal probabilities. Niemi et al j2H] assume n to be an odd 
number and Garman et al jU] indicate the relation in a footnote. However, 
they do not present a rigorous statistical proof that the limit is exactly a mul- 
tivariate normal probability for any general culture probabilities. They state 
that the multivariate normality is achieved by the fact that the multinomial 
distribution tends to the multivariate normal distribution. The multivariate 
normality however, is the consequence of the fact that a set of linear com- 
binations of Ni values, defined by equation (J2J), tends to the multivariate 
normal distribution by the application of multivariate central limit theorem, 
as we proceed to show here. In this section we show a rigorous statistical 
derivation of the limit which covers the cases when n is odd or even. 

Most researchers in this area have focused on the odd n (number of voters) 
case. This is because ties do not occur while determining the majority rule; 
i.e., for odd number of voters n, the number of voters preferring Ci to Cj 
will never equal the voters preferring Cj to Ci. However, if n is even there 
is a positive probability of having a tie. We show, in what follows, that the 
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probability of a tie for n even tends to zero as n — > oo. Let p^j = /]p. 



r 

r 



the with summation extending over the preference rankings in which C{ is 
preferred over Cj and let T denote the total number of voters preferring d 
to Cj. Then, it is well known from Multinomial distribution theory that T 
is a Binomial variable with n and as parameters. The probability of a tie 
in the majority determination between Ci and Cj is 

(n\ n . , n 



,2/ 

n! 



:[P^(1-Pii)] 2 (4) 



We know by Stirling's approximation Page 130, that 
The symbol = means " ! — i — > 1 as n — > oo. Applying this to equation 

\/27re~ n n n 2 



(0) and noting that — Pij) < | for all < p^ < 1, we can verify that 

„ i on+l 



< 



1 2 



71+1 



V2mr 4z 
^2 1 



— > as n — > oo 

7T 

For a given voter v, for i = 1, • • • , m, j= 1, • • • , m, % ^ j, define 

X(i,j,v) = 1 if CiPCj 

= -1 if C j PC l (5) 

n 

For the voter profile (N x , N 2 , ■ ■ - , N K ), QMCj if £x(z,j» > 1 . The 

«=i 

above representation has been used in Bell (0) and Gehrlein (|12j). Note 
that 

n 

v) = a UiJ) N l + a 2 ,(i,j)N 2 H h a K> ^N K 

v=l 
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We have, 

Xij = E(X(iJ,v)) 

= ai,(i,j)Pi + a 2 ,(i,j)P2 + ■■■ + a K ,(i,j)PK (6) 
Var(X(iJ,v)) = 1 - X% (7) 

We also need to obtain the correlation coefficient between X(i,j,v) and 
X(i, l,v), i ^ j, i ^ I and j ^ I. To express correlation, we first define 

a r,(ij),i — 1 if in the preference ranking corresponding to iV r , r = 1, • • • , K is 
CiPCj and QPCi 
or CjPd and QPd 
= —1 Otherwise (8) 

Covariance is defined by the following equation. 

Covariance(X(i,j,v),X(i,l,v)) = E[X(i, j,v)X(i,l,v)] — XijXu 

The correlation between X(i,j,v) and X(i,l,v) is given by 

(i) E(X(i,j,v).X(i,l,v)) - XjjXg 
K •> = = j = l,---,m, 

I — 1, • • • , to and j ^ I, j ^ i, I ^ i. 

K 

^ a r,(i,j),lPr — XijXu 
r=l 



For a given i = 1, 2, • • • , to, let Ri denote the (m—l)x(m — 1) correlation 
f). Note that R$ 



matrix (R^i). Note that R^ = l,j = 1, 2, • • • , to; j ^ i. 



Let 



V 1 - 

and Z„(^) = 

= ±^^M,J = l,-,m,iH (9) 

„=1 y/nJl — Xij 



Note that the (m — 1) random variables y(i,j,v),j = 1, • • • ,m, j ^ i have 
a joint distribution with zero mean vector, unit variances and Correlation 
matrix R4. 

From (jHJ), the probability that CiMCj is given by 



1 



y/nXi 



\v=l 



1 - Xfj 



(10) 



The probability that Ci is the Condorcet winner is therefore given by the 
joint probability 



P 



Z n (l,j) > ~ 



1-AJ- 



(11) 



By the multivariate central limit theorem (See Cramer [Oj, Page 112) 
applied to independent summands Z n (i,j),j = l,---,m, j ^ i, the joint 
distribution of 

{Z n (i, 1), ■ • • , Z n (i, i - 1), Z n (i, i + !),■■■, Z n (i, m)) 

tends as n — > oo to the (m — 1) dimensional multivariate normal distribution 
with zero mean vector, unit variances and correlation matrix Ri 
Let us write 

-oo if Xij > 
oo if \^ < 
if X^ = 



(if)- 



Denote by L(hi, hi, • • • , h m -i\ R) the multivariate normal probability 

P{Z{ > h U Z* > h 2 , ■ ■ ■ , Z* m _ x > hm-1] R) 

where (Z*, Z%, • • • , ^-i) ^as the multivariate normal distribution with zero 
mean vector, unit variances and Correlation matrix R. Then we have 



lim P 

n— >oo 



Z " > 'l-AJ = l ' 



L(Sn, Si2, • • • , 5t,i-lj ^i,i+l> ' ' ' 5 ^imj -Ri)- 



(12) 
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The limiting probability of a Condorcet winner is then 

m 

P(oo,m) = ^2L(8 il ,5i2,- " ,8^-1,5^+!,' ,8 im ;Ri). (13) 

i=l 

Equation (JTHj) shows that the limiting Condorcet winner probability is 
obtained as the sum of m L-functions. The arguments of the L functions are 
either +00 or —00 or 0. If at least one of the arguments in a L function is +00, 
the function value is by the definition of the L-function. We can also drop 
all the arguments with value —00 and calculate the reduced L-function from 
the marginal multivariate normal distribution. It is therefore interesting to 
note that even for arbitrary probability distribution p^s, at worst we have to 
evaluate positive orthant probabilities of a multivariate normal distribution. 
We shall give examples for the case m = 3 and m = 4 in later sections and 
provide a scheme to obtain the exact limiting probability expressions for a 
Condorcet winner for any choice of p^s. 

Equation (|13|) can also be specialized to the Dual Culture(DC) Model. A 
preference rank order is the dual of another ranking if each can be obtained 
by reversing the order of the other. For example, if m — 4, rank order 
C1C4C3C2 is dual to C2C3C4C1. The DC model assigns Pi such that Pi = pj 
if i is dual to j. In this case it is clear from the definition of that all the 
Xij = 0. This, in turn, implies that 8ij = 0. Equation (|13|) shows then that 
the limiting probability of a Condorcet winner under DC models 

m 

$X0,0,0,. -.,0;^) 

1=1 

Note that in each L function, there are (m — 1) zeroes. Gehrlein jTT] gives 
the exact expressions for m = 3 and m = 4. 

3 Condorcet Winner: Some Special Cases 

First we assume that all preference rankings for a voter are equally likely. 
This implies that Pi = 1/K, for i = 1,---,K. For this case, symmetry 
conditions imply that = for i — 1, • • • , m , j — 1, • - • , m and i 7^ j 
=^> 5ij = 0. For this case, it can be verified from ^3] that Ri = R* where 
R* is the (m — 1) x (m — 1) equi correlated matrix with equal correlation value 
= 1/3. The limiting probability of a Condorcet winner from (jl3|) is given by 

P(oo,to) = mL(0, 0, • • • , 0; R*) (14) 
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This expression for m = 3 was obtained by Gilbaud ([Ej) as §+§^ arcsin(~) 
Garman and Kamien [0] gave the expression for m = 4 without giving details. 
For m > 5, several researchers have obtained approximate values for expres- 
sion ([Hj) based on simulation studies or by quadrature methods. See for 
example Gehrlein J3], Jones et al and Niemi et al j2H]- The expression 



(|14jl shows that the limiting probability for m candidates can be obtained 
by an evaluation of positive orthant probabilities for a (m — 1) dimensional 
multivariate normal distribution with zero mean vector, unit variances and 
equicorrelated correlation matrix with all correlations equal to 1/3. Several 
papers in the statistical literature deal with such evaluations. See, for exam- 
ple, David and Mallows |7j, Sondhi[26j, Placket [21], Rub en [25]. Johnson and 
Kotz [IB] and Bacon |XJ . Bacon [T] derives explicit expressions for the positive 
orthant probabilities for a few small values of m for a multivariate normal 
distribution with equicorrelated correlation matrix with all correlations equal 
to p. He also derives a recursive relation to calculate the probabilities for 
successive values of m. Except for a few small values of m, such expressions 
involve the evaluation of multiple integrals without known closed form ex- 
pressions. Denote the expression (fT4*j) by P(m). Using the expressions from 
Bacon pp, Gehrlein ^U] obtains the values of P(m) for some small values of 



m as 



m = 3, P(3) = f + £ arcsin(f ; 



m 



4, P(4) = f[l + farcsin(i)] 



m = 5, P(5) = ^[1 + ±2 arcsin(i) + ^ g arcsin' A 1 ,/A 



m 



6, P(6) = ^[1 + f arcsin(|) + J$ arcsin( A ^ dA 



■ l+2A^ 



m = 7, P(7) = Ul + f arcsin(±) + arcsin' A ^ "' A 



64L 1 7T ^3; i n 2 ju " u - VA ""Vi+2A' 
720 rf ri+fe a^c^r A \ dA d/i 1 



+ ^^/o 1+2M arcsin( r+2A , 

Note however, that some of these expressions involve integrals. 

Denote L(0, 0, ■ • • , 0, p) by L m _i(p) when there are m — 1 zeros in the 
argument of L and p is the common correlation. From the results of Sampford 
(See Moran [21]) we have 

1 r°° ■> 1 
L m ^{p) = — e-* (l-$(at)) m -Vt 

where $(.) is the cumulative distribution function of the standard normal 
distribution and a = For the IC model, p = | and hence a = 1. Thus 

1— p ' " 3 
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J m— 1 



(0. 



6 WIT J-oo 



(15) 



It is also known that L m _i(p) for m being even can be obtained recursively 
from Li(p) for i < m — 2: 



Lm-i(p) = ^ 



— m 



1) 1 m-l\ _ . , 
— + 1+ 2 U 2 (p)- 



r 



^ 3 (p) + • • • + 



1 



ym — 2, 
(16) 

See Johnson and Kotz ^B] Pages 48-52 for more details. We obtained the 
expression for m = 6 using the recursion formula. 

For m > 7, Bacon pQ provides an approximation believed to be reasonably 
close to the true value. Figure [T] shows the plot of P(m) for increasing values 
of m based on these approximations. May jTH] has shown that P(m) — > 0, as 
m — > oo at the rate of — . His proof is based on saddle point approximation 
of an integral and it appears to be complicated and difficult to understand. 
We present a direct proof of this result. 

From (fTl^) and (fTB^l . we have 



\L 



m-2 



(p) 



P(m) 



m 

V2m 



e-* (1 - $(t)) m_1 dt 

-co 



do 
oo 



0(t)(l-$(t)) m_1 e-5rdt 



where 



P(m) 



V2: 



m 



m 



1 s ^£ 

e 2 



m j-oo 



[1 - $(t)] m te"2-rft 



-V2 /°° [1 - $(t)] m e^te^cft 

J —CO 



Applying Cauchy-Schwarz's inequality 



< 



f\t)dt\ 



f(t)g(t)dt 
to the above equation, we get 

Pirn) < V2{ f°° (1 - ®(t)) 2m e^dt 



g 2 (t)dtj 



(17) 



Ve—dt 
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Note that 

^ (1 - $ (t)) 2m e^dt 

J —oo 



/OO _ t 2 
t 2 e—dt 
-oo 

Thus 

P{m) < 2ttV2 

(2m +1)3 

We proved that the limiting probability Pirn) of a Condor cet winner 
tends to zero as m — > oo. We show in addition thatP(m) is monotonically 
decreasing in m. This settles Kelly's ^H] and Buckley and Westen's j3] 
conjecture for large values of n. This monotonicity property, we believe, is 
new. 

Proposition: P{m + 1) < P{m) 

Proof: From Equations ()14j) and (|15jl . we have to show that 



= V2^ r (1 - $(t)) 2m 0(t)dt 

J— oo 

~~ 2m + 1 

roo 

= V2n / t 2 0(t)dt 

J — oo 

1 



(ra + 1) 



7T 



/°° e" t2 (l - $(t)) m ^ <-^=r e-* 2 (l - ^(t))" 1 - 1 ^ (18) 

J —oo \l TT J —oo 



( ffl -I— 1 ^ /*oo „ /-OO .2 

1 ' ' e" 4 (1 - §(t)) m dt = V2(m + 1) / e^-(l-$(t))">(t)cft 



'7T J-oo 

(19) 

where <f>(t) = pdf of standard normal distribution. The right side of equation 
(fT9j) equals after integrating by parts 

-V2 [°° te^r(l-$(t)) m+1 dt 

J —oo 

Similarly the right side of equation (fTHJ) equals 



-V2 [°° te^(l-$(t)) m dt 

J —oo 
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Thus it is enough to show 

te 4- ( (i _ $(t)) m+1 - (i - $(t)) m ) dt > o 



Split the region of integration into /"^ + J °°,put t = —y in the first integral 
and after simplification and noting 1 — $(— y) = $(?/), it reduces to show 
that 



(1 - <%)) m+1 - (1 - Hy)) m ) - ($ m+1 (y) - $ m (y))] ^(»)dy > 

(20) 

We have written $ m+1 (y) for $(?/) m+1 and $ m (y) for $(y) m We will 
show that 

($m (y) _ (1 _ $(y)r) _ ^™+i (y) _ (1 _ > o for y > (21) 

Write $(y) = z and note that the left hand side of equation (}2"T]) equals 



(2z-r 



z m ' x + z m ' 2 (l -z) + ■■■ + (!- z) m - 1 



-{2z-l) 



z m + z m - 1 (l-z) + 



■■ + U-z) 



(22) 



Note that z m ~ x - z m_1 (l - z) = z m . Further 



Observe that 2z — 1 > for y > 0. 

)=*• 

™~ 2 (1 - z ) 2 = Z ™-\1 - z ) 



z m - 2 (l-z) 



z 



and so on. 

Hence the expression in the square brackets of equation (|20p 

= z m -\l -z) + z m ~ 2 (l -z) 2 + --- + z(l- z) m ~ l > 

Thus the integrand in equation (}2T)|) is >0 and hence the left side of equation 
(HD is > 0. 



4 Condorcet Winner: Special Case m = 3 
candidates 

Niemi et al (2H| calculated the limiting probability of a Condorcet winner 
values for n = 1..49, utilizing the table of Ruben [23]. Ruben used quadra- 
ture methods and recursive formulas to obtain these values. May JH] and 
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Garman et al discussed the case of 3 candidates to obtain the limiting 
probability of a Condorcet winner. Gehrlein ^T] discusses several possibili- 
ties that may arise in obtaining the value for m = 3. We give a compact and 
complete Table for all possible p,'s utilizing the results of Section 2 and (|T3*j) 
of this paper. Suppose the 6 rankings are 



Rankings 


Number of Voters 


Probabilities 


C1C2C3 


iVi 


Pi 


C1C3C2 


iV 2 


P2 


C2C1C3 


N 3 


P3 


C2C3C1 


N 4 


Pa 


C3C1C2 


N 5 


P5 


C3C2C1 


N 6 


P6 
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It can be verified that 



Al2 


= P1+P2+P5 


~P3 


-Pi 


-Ps 


Al3 


= P1+P2+P3 


-Pa 


-P5 


-P6 


A21 


= P3+P4+P6 


-Pi 


-P2 


~Pb 


^23 


= P1+P3+P4 


-P2 


~P5 


-Ps 


A31 


= P4+P5+P6 


-Pi 


-P2 


~P3 


A32 


= P2+P5+P6 


-Pi 


~P3 


-Pa 



The three correlation matrices Ri, R2, R3 are given by: 



Ri 



( 1 

(Pl +P2 +P4 +P6 ~P3 ~Pb ) ~ A 1 2 A 13 

V 



(Pl+P2+P4+P6-P3-Ps)-Al2Ai3 

V /(1-A? 2 )(l-A? 3 ) 



V /(1-Af 2 )(l-Af 3 ) 

The (1,2) entry of Ri is R\ z . 

( 1 

(P2+P3+P4+P5-Pl-P(i)-A2lA23 

V 



Ri 



v /(l-Ai 1 )(l-Al 3 ) 

The (1,2) entry of R 2 is R 2 n . 



R3 



1 



(Pl +P3 +P5 +P6 -P2-P4)-A3lA32 



(P2 +P3 +P4 +P5 -P!-P6)"A2lA23 

v /(i-a| 1 )(i-a2 3 ) 

1 



(Pl +P3 +P5 +PG -P2-P4)-A3lA 32 

V /(1-A| 1 )(l-A§ 2 ) 
1 



V /(1-A| 1 )(l-A| 2 ) 

The (1,2) entry of -R3 is i?f 2 - The probability of a Condorcet winner is 

£(<*i2, *i 3 , #23) + L&i, S23, R 2 13 ) + Ufa, $32, R 3 u)- (23) 



By considering the 27 possibilities with which A12, A13 and A23 can take zero, 
positive and negative values, we can simplify the expression further. Table 
[0 shows the 27 possible cases and the corresponding limiting value of the 
probability of a Condorcet winner. A set of different values is obtained for 
varying values of p^s. It is also to be noted that for a few alternatives the 
values of and 1 are obtained for the limiting probability of a Condorcet 
winner. See alternatives 7 and 17 in Tabled Table ^ provides the reference 
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table to calculate the exact limiting probability for any choice of Pi,i = 
1, • • • , 6. One of the 27 possibilities will apply in a given situation and one 
can read off or calculate the value. It is seen from Tabled that only in two 
cases out of 27, the limiting probability is definitely 0. In 12 cases (more 
than 40%) the probability is 1 and in almost all cases the probability will be 
at least 50%. These observations are consistent with often stated conclusion 
that for small number of candidates and large number of votes, the Condercet 
scheme will produce a winner. 

It can be verified that case 7, for example, is reached when pi = P2 = 
P5 = §; P3 = | and p 6 = jq. Case 17 is reached when p x = ^,p 2 = ^,^3 = 
iI>P4 = if) Ps = 22 anc ^ ^6 = 0. For another set of examples for m = 3, see 
May [E!, Garman et al and Meyer et al |2Uj. 

For example, the trivial set of probabilities Pi = l,pj = 0,j ^ i will 
produce a Condorcet winner for every n and hence the limit is also 1 as n — > 
oo. Choose the probabilities such that C1C2C3C4 ■ ■ ■ C m , C^C^CiC^ ■ ■ ■ C m , 
C3C1C2C4, ■ ■ ■ C m have probabilities | each and the rest of the rankings have 
probabilities 0. It is easily verified that A i2 = |, A i3 = A23 = |, Xij = 
X 2 j = A 3 j = 1 for all j > 4. This implies that <5i 3 = 00, 612 = 623 = —00, 5ij = 
$2j = 833 = —00, for j > 4. From equation (|T2"j). it follows immediately that 
in every L term, one of the arguments is a +00. This implies that the limiting 
Condorcet winner probability is 0. Thus for every m, we can choose Pi so 
that the limiting Condorcet winner is 0. 

Gillettjni] conjectured that the probability of a Condorcet winner P(3, n; p) 
for m = 3 candidates and any n, the number of voters, is minimized for 
Pi = Pa = P5 = l,P2 = P3 = P& = or p x = p 4 = p 5 = 0,p 2 = p 3 = Pe = |- 
Buckley proved this conjecture for the special case n = 3. We will show 
the validity of the conjecture for general n. 

The probability of candidate C\ being the majority winner is 

Ni + N 2 + N 3 > - and N x + A^ 2 + iV 5 > -J 

= P^N l + N 2 >^- mtn(N 3l iV 5 ) j 

>pl [ N 1 + N 2 >^Y 

since A^3 and A^ 4 are non-negative random variables. The above expression 
is greater than or equal to 
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P jiVi > §}, again for the same reason with regard to N 2 . N\ is dis- 
tributed as a Binomial variable with n and p\. 

p(N l > T ^=l-B{k ] n, Pl ) 

where B(k;n,pi) is the Binomial cumulative distribution with parameters n 
and pi. k = if n is odd and k — | if n is even. From the relation of 
B(k;n,pi) with Incomplete Beta function, see Feller(1957) page 127, jS], 

P( Nl >^)=n( K n ~ k 1 "j jT t fc (l - tr^dt (24) 

Continuing similar arguments for C 2 and C3 being majority winners, we 
have 



P(3, n;p) > n 



k 



ppi fP4 rP5 

/ t k {l-t) n - k - 1 dt+ / t fc (l - t) n_fc-1 dt + / t k (l-t) n - k - l dt 
Jo Jo Jo 



Consider minimizing the following function 

tt(Pl,P4,P6) = [ P1 t k (l-t) n - k - 1 dt+ fV(l-t) 

Jo Jo 

+ / t k a-t) 

Jo 



n - k - x dt 



>[ dl 



subject to p\ + P4 + p 5 = c 
Pi > 0,p 4 > 0,p 5 > where c = 1 - (p 2 + P3 + P&) (25) 



We proceed to show that an optimal solution to the problem defined by 
(|2~5Jl is given by p\ = p 4 = p 5 = |. 

it(pi,P4,p 5 ) is a separable function in p l5 p 4 and p 5 . It can be verified 
that 

9 2 n( Pl ,p 4 ,p 5 ) = _ Pi n-,- 2) (Jfc _ (n _ fi = 1;4;5 (26) 



2 



dpi 

k 



For ^ < -^r, i = 1,4, 5, -f^| > =>• u(pi,p4,p 5 ) is convex in pj in the region 



of p^s. 

In order that p^s satisfy the constraints of equation (j23j) . we cannot have 
two of the Pi's are > -^-r with at least one of them > since their sum 

r 1 — n— 1 n— 1' 

will be greater than > 1 for both n even and n odd. If for any feasible 
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solution, a single Pi > -^r, which implies that the other two < -^j, the 
Hessian matrix will have its determinant a negative value because of equation 
([26)1 . Note that the Hessian is a diagonal matrix. Thus such a solution cannot 
be a local minimum. Therefore we can confine our solution to the region in 
which all the p^s < It is easy to verify that p\ = p 4 = p 5 = | satisfies 
the Kuhn- Tucker condition for the minimum. Note also that § < , fc , , since 

3 (n— 1) ' 

> 1. Since ^(^1,^4,^5) is convex in this region, the solution is a global 
minimum of u(pi,p4,p 5 ). 

Fix now pi = p 4 = p 5 = | and repeat the whole argument with P2,P3, Pq as 
variables and consider the minimization problem of minimizing the following 
function. 

w(p 2 ,p 3}P6 ) = / t k (l-t) n - k - 1 dt+ t\l-t) M dt 
Jo Jo 

+ r t k {\-t) n ~ k - i dt 

Jo 

subject to pi + pz + Pq = 1 — c 
P2 > 0,P3 > 0, p e > where c = (pi + p 4 + Ps) (27) 

As before, the minimum solution of equation (J27|) is pi = p% = p$ = ^ 
for pi — p 4 = p 5 — |. From equation (j2ZJ), we see that the least value 
of w(pz,P3,pq) is when c = 1. Thus the minimum solution is pi = p± = 
P5 = |; £>2 = £>3 = P6 = 0. Note that our arguments show that p% = ps = 
P6 — §j Pi — — £>5 — is also optimal. Denote by p* the solution 
Pi* = p^ = ps* = §; P2* = Pz* = Pa* = 0. Let u* (pi* , p 4 * , p 5 *) be the value 
of the objective function of (J25j) . We have shown that for any choice of p 

P(3,n;p) > fc Ju*(px*,p4*,p 5 *) 

= 3(1- B(k;n,h), from fEJ 

We will now show that P(3,n; p*) = 3(1 — B(k; n, |)) which concludes 
that p* yields the minimum value of P(3,n;p). Since p 2 * = Pz* = p$* = 0, 
we have 

71 11 

P(3,n;p*) = P( Nl >-,N 1 + N 5 >-) 

Tl Tl 

+ P(N 4 > -N i + N 1 > -) 
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+ P(N 5 >^,N 5 + N A >^) 

= P(Ni>^) + P(N4>l) + P(N 5 >l) 

= 3(1 -B(k;n,^), 

since Ni,N A ,N 5 are identically Binomial distributed with parameters n and 
1 

3' 

This establishes the conjecture. Note that the minimum value for the 
probability of a Condorcet winner is given by 



or an equivalent integral representation 

'n — 1 



3n\ 



k 



Ht k (i-t) n -"- 1 dt 

Jo 



The former is easier to evaluate for small values of n and employ a normal 
approximation to the Binomial for large n. For example, when n — 3, the 
minimum probability is 

which is about 77%. For n — 4, the minimum probability of a Condorcet 
winner is 

*((i)a)'(V\ + '* n ^ 1 



3/ \SJ \3J \4J\3J J 3 

which drops down to 33%. 

The foregoing arguments can be generalized to m > 3 candidates. Con- 
sider the probability distribution p* on to! preference rank orders given by 
assigning probability of ^ to each of the cyclic rank orders 

rank order 1, 2, • • • , (m — 1), to probability pi 
rank order 2, 3, • • • , to, 1 probability P2 
rank order 3, 4, • • • , 1, 2 probability p 3 

rank order to, 1, • • • , (to — 2), (to — 1) probability p m (28) 
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We will establish that this solution yields the minimum for the probability of 
a Condorcet winner P(m, n; p) for a general probability distribution p. We 
use the same approach that was used for m = 3 case. It is shown first along 
similar lines that 



77 Tt Tt 

P(m,n;p) > P(N 1 >-) + P(N 2 >-) + ■■■ + P(N rn > -) 



n 



'n-V 
k 



i=l 



E/ t\i-t) 



n—k—l 



(It 



From an argument similar to the case m = 3, it is shown that 

P(m, n; p) > mil — B(k;n, — )). 

m 



(29) 
(30) 

(31) 



Next it follows that from (f2*Kj) and the definition of a Condorcet winner, 
P(m, n; p*) for p* given by pi* = — , for i = 1, ■ ■ • , m and pi* = for 
i > (m + 1) the probability of C\ winning is 

77 77 

P ( Nl + N 3 + . . . + Nm > Nl + iv 4 + ■ ■ ■ + N m > -, 

-^ 1 >-) 



P(iVi> 



n , 



Similarly it is verified that the probability of C 2 winning is 



n 



P(N t + N 2 + N 4 + ■ ■ ■ + N m > -, Nt + N 2 + N 5 + ■ ■ ■ + N m > 

■■-,N 2 > 



and so on. Hence 
P(m, 77; p*) 



n 



n , 



n 

2' 

n , 



n , 



P(N, >-) + P(N 2 >-) + ... + P(N m > - 



n , 



P(N 2 > -) 



m(l — B{k] 77, — )) 



777 



which equals the bound in (|31|). Hence p* yields the minimum for the prob- 
ability of a Condorcet winner. Table El gives the minimum probability of a 
Condorcet winner for selected values of m and 77. 
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5 Condorcet Winner: Special Case m=4 can- 
didates 



For this case, we show that it is possible to obtain, for any choice of i = 
1, 2, • • ■ , 24, exact limiting probability of a Condorcet winner. We first cal- 
culate the six quantities Ai 2 , Ai 3 , A 14, A23, A 2 4 and A 34 from (jUJ). Then define 
812, 813, £14, ^23, ^24 and £34 as stated in Section 2. Also calculate the correla- 
tion matrices Ri,i = 1,2,3,4. In this case, there are 3 6 = 729 possibilities 
for the Ay's to be > or < or = 0. It will be unwieldy to prepare a Table 
containing all these 729 cases. We show how the limiting probability can 
be calculated for any given situation. The limiting probability of Condorcet 
winner follows from (|13|). 

P(oo, 4) = L(5 12 , 5i3, 6 U ; Ri) + L(S 2 i, 5 23 , £24; #2) 
+L(5 3 i, 5 32 , £34; R3) + £(<$41, &t2, (^43; R4) 

Consider the calculation of L(5\ 2 , £13, 8u\ R\). Calculations for other L- 
factors are similar. To obtain this, we need at the worst to calculate L(0, 0, 0; R{) 
or L(0, 0, Rjj) or L(0). Recall the definition of R\, as appropriate i, j element 
of R\. Note that L(0) = $(0) = ~, where $ is the cdf of the standard normal 
distribution. Bacon 1] gives exact expressions for those L functions as 

L(0,0,4) = - A + 1-arcstniR^) 



L(0,0,0,R 



2 

1 H — (arcsin(Rl 3 ) + arcsin(Rl 4 ) + arcsin(R\ A 



The calculation of this L for other possible arguments is by inspection or it 
will reduce to one of the above 3 cases with arguments equal to zero. 

This example illustrates the comment made immediately after equation 
(fT3*|) that the only calculations, if any, involve calculating the positive orthant 
probabilities of appropriate multivariate normal distribution. Thus for any 
one of the 729 possibilities that is reached in a given case, we can calculate 
the exact limiting probability of a Condorcet winner. 

6 Conclusion 

In this paper, we provide a rigorous mathematical treatment for calculating 
the limiting probability of a Condorcet winner for any number of candidates 
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and with arbitrary rank order probabilities, when the voters are independent. 
We show further that the limiting probability depends only on the positive 
orthant probabilities of an appropriate multinomial distribution even for ar- 
bitrary probabilities of the preference rank order probabilities. We provide 
approximate limiting probabilities for a Condorcet winner for any number of 
candidates and equal rank order probabilities. We provide a compact and 
complete Table for the limiting probability of a Condorcet winner with 3 
candidates and arbitrary rank order probabilities. We present a simple proof 
of a result of May [19] to show that the limiting probability of a Condorcet 
winner tends to zero as the number of candidates tends to infinity. We also 
present a scheme for calculating the exact limiting probability for m = 4 
candidates and arbitrary probability distributions on the preference rank or- 
ders. We show for the first time that the limiting probability of a Condorcet 
winner for any given number of candidates m is monotone decreasing in m 
for the equally likely case. This, in turn, settles the conjectures of Kelly [TH] 
and Buckley and Westen [3] for the case n —>■ oo. We prove the validity of 
Gillett's jin] conjecture on the minimum value of the probability of a Con- 
dorcet winner for m = 3 and any n. Extending the result we obtain the 
minimum value and the minimum solution for any n and m. 
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Number 


Al2 


A13 


A23 


Probability of a Winner 


1 











\ + ^ axcsin(i^ 3 ) + \ + ^ arcsin^) + | + ^ arcsin( J R? 2 ) 


2 








+ 


k + J + h arcsin(i?y 


3 








- 


\ + J + h arcsin(i^) 


4 





+ 





1 + 1 + 2^ arcsin(^ 3 ) 


5 
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1 + 1 + 2^ arcsin(i?y 


6 





+ 


- 


1 
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7 





+ 
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1 


8 





- 


+ 


1 
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9 





- 


- 


1 


10 


+ 





+ 


1 

2 


11 


+ 


+ 





1 


12 


+ 





- 


1 


13 


+ 


- 





1 

2 


14 


+ 
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1 


15 


+ 
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16 
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+ 
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17 


+ 




+ 





18 


+ 








2" + I + 2^ arcsin( J Rf 2 ) 


19 










5 + j + £ arcsin(i?? 2 ) 


20 
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21 









1 
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22 
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23 
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24 




+ 







25 
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26 
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Table 1: Limiting Condorcet Winner Probabilities for m = 3 candidates 
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n 


m = 3 


m = 4 


m = 5 


m — 10 


3 


0.7778 


0.6250 


0.5200 


0.2800 


4 


0.3333 


0.2031 


0.1360 


0.0370 


5 


0.6296 


0.4141 


0.2896 


0.0856 


6 


0.3004 


0.1504 


0.0848 


0.0127 


7 


0.5199 


0.2822 


0.1667 


0.0273 


8 


0.2638 


0.1092 


0.0520 


0.0043 


9 


0.4345 


0.1957 


0.0979 


0.0089 


10 


0.2297 


0.0789 


0.0318 


0.0015 


19 


0.1943 


0.0356 


0.0079 


0.0000 


20 


0.1129 


0.0158 


0.0028 


0.0000 


29 


0.0934 


0.0071 


0.0007 


0.0000 


30 


0.0564 


0.0033 


0.0003 


0.0000 


50 


0.0148 


0.0002 


0.0000 


0.0000 


51 


0.0207 


0.0002 


0.0000 


0.0000 


100 


0.0006 


0.0000 


0.0000 


0.0000 


101 


0.0008 


0.0000 


0.0000 


0.0000 



Table 2: Minimum Condorcet Winner Probabilities for some values of m and 

n 
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